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Introduction 


A  review  of  the  basic  derivations  leading  to  the  standard  electromagnetic  wave  equations 
reveals  a  number  of  errors.  These  errors  and  the  revisions  required  to  correct  the  basic 
formalism  are  the  subject  of  this  paper.  The  present  section  lists  the  basic  Maxwell’s 
equations  and  provides  an  outline  of  the  standard  wave  equation  derivations  for  later 
reference.  Subscripts  and  superscripts  are  used  in  an  attempt  to  provide  precision  to  the 
definition  of  the  key  variables.  Additional  justification  for  some  of  the  notation  used  here 
will  be  provided  in  the  course  of  the  paper. 


V  •  E*  =  p/s . 

(1) 

V  x  Ec  =  0  . 

(2) 

V  •  B  =  0  . 

(3) 

VxEj  =  -fiB/fit . 

(4) 

VxB  =  pJT  -l-psfiEc/fit  +  psfiEj/fit . 

(5) 

s  is  the  electrical  permittivity  and  p  is  the  magnetic  permeability. 
Eq.  (1)  is  Gauss’  law  for  a  static,  or  quasi-static  (field  fluctuations  are  transmitted 
instantaneously)  Coulomb  field  E®  generated  by  a  Coulomb  charge  distribution  p.  It  is  a 

consequence  of  the  inverse  square  law  of  the  Coulomb  field  of  a  point  charge.  Eq.  (2) 
reflects  the  fact  that  the  general  Coulomb  fields,  including  static  and  dynamic  fields,  are 
conservative  and  can  be  expressed  as  the  gradient  of  a  scalar  Coulomb  potential, 
Ec  =-Vcpc.  (In  the  foregoing,  Ec  generally  refers  to  dynamic  fields.)  Equation  (3) 

states  that  B  is  solenoidal.  Equation  (4)  is  Faraday’s  law.  Given  the  magnetic  vector 
potential,  A,  where  B  =  V  x  A ,  it  follows  from  Eq.  (4)  and  the  assumption  that  A  is  the 
sole  source  of  Et ,  thatEj  =  -  fiA Id  t . 

Eq.  (5)  is  obtained  from  Ampere’s  circuital  law.  The  right  hand  side  is  the  sum  of  true 
currents,  JT  and  the  displacement  current,  JD .  JD  is  generally  the  sum  of  contributions 

from  time  derivatives  of  Ec  and  E, .  Self-consistency  requires  a  solenoidal  net  current 
source  (sum  of  JT  and  JD)  for  B. 

In  the  standard  derivation  for  the  wave  equations  in  terms  of  the  potentials,  one  uses  A 
and  (pc  in  Eq.  (5),  to  obtain, 


V x V x A  =  V(V •  A) - V2A  =  pJT  -psfiV(pc/fit-psfi2A/fit2  (6) 

Rearranging  terms, 

V(V  •  A  +  psfi(pc  /fit)  =  V2A  -psS2A/fit2  +  pJT .  (7) 
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The  two  terms  in  parentheses  on  the  left  hand  side  of  Eq.  (7)  are  assumed  to  be 
independent  and  unrelated,  with  V  •  A  treated  as  indeterminate.  The  rationale  is  that, 
sinceB  =  VxA,  any  gradient  function  can  be  added  to  A  without  affecting  B.  If 
somehow  the  left  hand  side  could  be  set  to  zero,  Eq.  (7)  would  represent  an 
inhomogeneous  wave  equation  with  source  tenn,  JT .  The  conventional  approach  for 

obtaining  a  wave  equation  from  Eq.  (7)  invokes  the  electromagnetic  gauge  which 
transforms  the  laboratory  system  of  variables  to  new  variables  so  that  the  left  hand  side  of 
the  transformed  version  of  Eq.  (7)  is  zero  (the  Lorenz  condition).  More  specifically,  the 
procedure  involves  application  of  the  gauge  transformation  function,  \\i ,  such  that 


A  -»  A’  +  Vvp 

(8) 

cpc  — >cpc  -  d\|/  /  <3t . 

(9) 

These  transformations  allow  arbitrary  alterations  to  the  variables  without  affecting  the 
total  electric  field,  E.  It  is  readily  seen  from  Eqs.  (8)  and  (9)  that 

E  =  Ec  +  E,  =  EC  +  Ej,  (10) 

where  Ec  =  -V(pc ,  E,  =  -  8AJD  t ,  Ec  =  -V(pc ,  and  Ej  =  -  <3A'/d  t . 

To  obtain  the  wave  equation,  first  substitute  Eqs.  (8)  and  (9)  into  the  left  hand  side  of  Eq. 
(7),  and  rearrange  to  give, 

V(V*A  +  ps5(pc /fit)  =  V(V  •  A  +  ps5(p/<9t  -  V2v|/  +  psd2v|//dt2)  .  (11) 


One  obtains  the  Lorenz  condition, 

V(V  •  A  +ps<9(pc /6t)=0,  (12) 

from  Eq.  (11),  by  selecting  the  function  vp  such  that  its  values  in  time  and  space  give  zero 
for  the  sum  of  terms  in  parentheses  on  right  hand  side  of  Eq.  (11),  i.e., 

v2^-d\/dt2  =  (V«A  +  S<pc/dt)  .  (13) 

If  we  now  return  to  Eq.  (7),  and  apply  Eqs.  (8)  and  (9),  and  the  Lorenz  condition  (Eq. 
(12)),  we  obtain, 


V2A'  -psa2A/fit2  +  pJT=0.  (14) 

Eq.  (14)  is  the  inhomogeneous  wave  equation  for  the  transfonned  vector  potential 
(primed  variables). 
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To  complete  the  task  of  solving  for  the  total  field  E,  one  needs  the  corresponding  scalar 
wave  equation  in  the  transfonned  variables.  This  requires  the  dynamic  form  of  Gauss’ 
law,  which  is  derived  from  the  total  electric  field, 

E  =  Ec+E;.  (15) 

The  key  step  here  is  to  assume  Maxwell’s  Eq.  (1)  holds  for  the  total  field,  E.  So, 
inserting  Eq.  (15)  into  Maxwell’s  Eq.  (1)  gives  the  familiar  result, 

V  •  (Ec  +  E'j)  =  p/s  .  (16) 

Note  that  Eq.  (16)  must  also  hold  for  the  unprimed  (untransformed)  variables  since  it 
involves  the  sum  of  the  two  fields. 

Expressed  in  terms  of  potentials,  Eq.  (16)  becomes, 

V  •  (-Vcpc  -  <9A’  /  fit)  =  p/e .  (17) 

Inserting  the  Lorenz  condition  (Eq.  (12))  into  Eq.  (17)  gives  the  scalar  wave  equation, 


V2(pc-psS?(pc/<3t2  =  -p/s.  (18) 

The  curl  and  time  derivative  of  the  solution  for  (14)  gives  B,  and  E,,  respectively.  The 
gradient  of  the  solution  for  Eq.  (18)  gives  Ec  .  Summing  Ec  and  E,  completes  the  task 
of  obtaining  the  total  electric  field  E  in  Eq.  (15),  which  is  independent  of  the  choice  of 
gauge  function. 

The  preceding  outlines  the  standard  formalism  for  the  elementary  derivation  of  the  wave 
equations.  The  remainder  of  this  paper  addresses  questions  regarding  the  validity  of  these 
standard  results  and  offers  alternatives  that  correct  a  variety  of  errors. 


The  Missing  Maxwell  Equation 

The  first  item  that  we  address  is  the  missing  Maxwell  equation  in  the  standard  formalism. 
Our  approach  maintains  the  distinctions  between  the  different  physical  variables  such  as 
E,  Ej ,  and  Ec  .  While  this  leads  to  a  proliferation  of  variables,  it  provides  more  clarity.  In 

the  literature,  the  same  symbol,  E,  is  used  for  all  electric  field  variables,  which  obscures 
the  physics  and  leads  to  errors,  as  we  will  illustrate  (see  also,  the  Appendix).  The  same 
problem  exists  with  the  various  components  of  the  vector  potential. 

The  fact  that  the  basic  set  of  Maxwell’s  equations  is  incomplete  is  demonstrated  by  an 
error  in  the  standard  derivation  of  the  wave  equation  for  the  electric  field.  The 
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conventional  approach  proceeds  as  follows  (see,  for  example,  Panofsky  and  Phillips, 
Chapter  11  [1],  Jackson,  Chapter  7  [2]  ,  Feynman  et  al,  Chapter  20  [3]  ):  take  the  curl  of 
both  sides  of  Eq.  (4),  which  gives, 

VxVxE,  =V(V*Ei)-V2Ei  =  -d(VxB)/St.  (19) 

From  Eq.  (19)  and  Maxwell’s  Eq.  (5)  we  have, 

V(V  •Ei)-V2Ei  =  pdJT/dt-psd2Ec/dt2- psd2E,/dt2.  (20) 

The  final  step  invokes  Maxwell’s  Eq.(l)  to  justify  setting  V  •  E ,  =  0  .  Treating  the  first 

two  terms  on  the  right  of  Eq.  (20)  as  external  source  terms  gives  inhomogeneous  wave 
equation, 

V2Ej-  ps^Ej  /fit2  =  pJT  +<92Ec  /fit2 .  (21) 

(We  derive  Eq.  (21)  later  from  the  vector  potential  wave  equation.)  Maxwell’s  Eq.  (1)  is 
usually  invoked  again  at  this  point  [1,  2,  and  3]  to  provide  the  standard  proof,  using 

V  •  Ej  =  0  ,  that  the  vector  E,  is  transverse  to  B  and  to  the  direction  of  wave  propagation. 

Glossing  over  distinctions  among  variables  is  the  likely  explanation  for  the  improper  use 
of  Maxwell’s  Eq.(l)  to  set  V*Ej  =0  .  The  variable  E,  refers  to  an  induced  field,  which 

is  obtained  from  the  vector  potential,  while  Maxwell’s  Eq.  (1)  holds  only  for  the  static 
Coulomb  field,  with  the  Coulomb  charge  density  as  its  source  term. 

It  follows  from  this  example  that  the  standard  wave  equation  for  E,  has  yet  to  be  derived 
properly.  One  needs  an  expression  for  V  •  E(  for  a  valid  derivation.  Another  direct 
indication  that  something  is  missing  from  the  basic  set  of  equations  is  that  the  dynamic 
form  of  Gauss’  law  (Eq.  (16))  is  supposed  to  give  an  expression  for  the  source  density 

V  •  E,  in  the  presence  of  a  dynamic  Coulomb  field,  while  no  such  expression  exists  for 
the  static,  quasi-static,  or  charge  free  cases. 

Given  that  the  validity  of  equation  (21)  is  not  in  doubt,  the  missing  Maxwell  equation 
must  be 


V  •  E,  =  0  .  (22) 

More  complete  justifications  for  Eq.  (22)  are  provided  in  the  following.  We  will  also 
show  that  Eq.  (22)  does  not  conflict  with  the  correct  form  of  the  dynamic  Gauss’  law. 
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Vector  Potential 


The  general  expression  for  the  vector  potential,  AG ,  due  to  a  general  current  density,  JG  , 
is  derived,  for  example,  in  Panofsky  and  Phillips  [1], 


Ag 


F  1 

r  JG(r') 

471  * 

'  r-r’ 

dv\ 


(23) 


A  key  point  here  is  that  AG  is  only  defined  within  an  arbitrary  function  that  has  a 

vanishing  curl.  As  shown  in  Panofsky  and  Phillips  [1],  any  three  dimensional  vector  is 
fully  characterized  by  giving  its  curl  and  its  divergence.  Generally,  if 

V  x  F  =  K  ,  (24) 


so  that  K  is  solenoidal,  and, 


V • F=s , 


(25) 


it  follows  that  F  may  be  expressed  as, 


F  =  -V  (pF  +  V  x  L , 


(26) 


where, 


and, 


9f 


1  I 

r  s(r')f 

An 

1-^ 

(27) 


L  = 


-LfK«dv. 

An J  r-r’ 


(28) 


Consider  a  general  vector  potential,  AG  where  V  x  AG  =  B .  An  expression  for 
V  •  AG(orcpA)  is  required  for  complete  characterization  of  A,  so  that,  generally, 


Ag  =-V(Pa  +VxLa, 


(29) 
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where 


L 


A 


1  fB(r') 

4  n  ■*  r-r’ 


dv'. 


(30) 


Similarly,  forEI?  we  have  VxE,  =  -SB/ft,  and,  allowing  the  possibility  of  a  source 
potential,  cpj ,  the  complete  E,  is  given  by, 


Ej  =  -Vtpj  +  VxLj,  (31) 

where,  invoking  Maxwell’s  Eq.  (4),  we  have 

1  4  n]  |r-rj 

Maxwell’s  Eq.  (4)  (the  original  Faraday’s  law)  gives  only  the  curl  ofEj.  The  divergence 
of  Ej  is  also  required  in  order  to  fully  characterize  this  variable.  In  other  words,  the 
variable  Et  is  undefined  so  long  as  V  •  E,  remains  undefined.  So,  we  have  returned  here 
to  the  case  of  the  missing  Maxwell  equation  (Eq.  (22)).  If  E;  were  arbitrary,  a  broad 
range  of  commonplace  physics  and  engineering  problems  involving  induced  fields  could 
not  be  addressed.  As  discussed  in  the  example  of  a  closed  circular  loop  in  the  Appendix, 
approximately  a  century  and  a  half  of  experience  with  applications  of  the  Faraday  law 
indicates  that  no  point  sources  exist  forEj.  Consequently,  the  required  equation  for  a 

closed  loop  of  current  is  V  •  E(  =  0  (Eq.  (22)).  This  equation  must  be  added  to  the  usual 
statement  of  Faraday’s  law  (Maxwell’s  Eq.  (4))  for  a  complete  expression  of  Faraday’s 
law. 

As  reviewed  in  the  Introduction,  the  standard  formalism  treats  V  •  AG  as  completely 
arbitrary.  In  view  of  the  above  discussion  of  the  Faraday  law,  this  approach  is  no  longer 
acceptable.  Consider  the  net  vector  potential  A.  Since  <9(V  •  A)/5  t=V  •  E(  =  0 ,  it  follows 
that  V  •  A  is  restricted  to  a  time  independent  function.  Consequently,  when  Eq.  (22) 
applies, 


8(V  •  A)/St  =  0.  (33) 

How  is  Eq.  (33)  reconciled  with  special  relativity  where  the  Lorenz  condition  (Eq.  (12)) 
is  generally  applied  to  both  A  and  A  ?  Again,  we  will  show  the  answer  lies  in  preserving 
distinctions  among  variables. 
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Alternative  Approach  to  the  Wave  Equations 


Retarded  Fields 


Among  the  key  equations  needed  for  our  purposes  are  the  retarded  Coulomb  and  vector 
potentials  of  a  moving  charge.  These  potentials  deal  with  the  essential  features  and  thus 
provide  a  more  substantive  basis  for  the  development  of  the  wave  equations  than  that 
used  in  the  standard  fonnalism.  A  further  advantage  is  that  these  equations  are  consistent 
with  the  requirements  of  special  relativity.  The  expressions  for  the  Coulomb  and  vector 
potentials  given  in  Feynman  et  al  [3]  for  a  charge  moving  along  the  x  axis  at  velocity  v 
are: 


(pc  = 


47isVf 


■V  /c" 


(x-vt)2 

l-v2/c‘ 


v  1/2 


(34) 


■+y2+z2 


and 


A.,  =■ 


(q/c2)  v 


47isVl-v2/c2 


(x-vt)2  2  2 

Vr+y+z 

l-v2/c2 


1/2 


A  =A_  =0. 


(35) 


J 


Recall  that  these  equations  reflect  the  fact  that  the  potentials  and  vector  fields  at  a  point  r 
are  given  by  the  moving  charge  and  current  located  at  a  point  r  (t  )  that  is  different  from 
the  present  source  position  r  (t)  because  of  the  finite  speed  of  light;  in  these  cases 
f  =  t-  r(t)-r  (t)  /c  .  The  term  -  ( v/c ) 2  originates  from  the  Lienard-Weichert  expression 

for  the  potential  and  reflects  the  apparent  elongation  of  the  charge  in  the  direction  of 
motion  due  to  retardation  effects  [3]. 


For  later  reference,  the  concept  of  retarded  fields  requires  a  modification  to  the  general 
expression  (Eq.  (23))  for  a  general  vector  potential  AG  due  to  a  general  current  source, 

JG(r  ,t  ) .  The  “retarded  solution”  version  [1,2]  of  Eq.  (23)  is: 


Pofr'^dv' 

r(t)-r(t) 


(36) 


where  the  source  term  in  square  brackets  is  evaluated  at  time  t  <  t . 
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The  Lorenz  Condition 


Our  approach  relies  upon  Eqs.  (34)  and  (35)  to  obtain  basic  relationships  among  the 
variables.  Eqs.  (34)  and  (35)  replace  the  vague  notion  of  “dynamic  fields”  with  precise 
definitions.  It  is  assumed  here  that  “dynamic  fields”  is  synonymous  with  retarded  fields. 
The  same  relationships  hold  for  any  charge  distributions,  in  view  of  the  principle  of 
superposition  for  sums  over  distributions  of  charges  q;  and  distributions  of  current 

elements,  q;V;  .  Note  that  Ax  =  vcpc  /  c2 . 

Eq.  (35)  is  the  vector  potential  associated  with  a  moving  charge  and  is  therefore 
associated  with  an  element  of  the  true  current,  JT .  We  now  label  this  true  current 

component  of  the  vector  potential  asAT,  where  generally  AT=ATXi+ATYj+ATZk  (  i,  j, 

and  k  are  the  unit  vectors).  There  are  further  contributions  to  the  net  vector  potential,  A, 
from  displacement  currents,  which  will  be  considered  when  we  discuss  the  vector  wave 
equation.  For  now,  we  examine  only  the  relationships  between  the  variables  pertaining  to 
true  currents.  The  first  item  is  that  the  relation  between  the  Coulomb  field  due  to  a 
moving  charge  and  the  vector  potential  due  to  the  current  due  to  that  moving  charge  is 
now  expressed  as  AT  =  vcpc/c2 . 


Next,  if  one  computes  the  divergence  of  Ax  ,  and  the  time  derivative  of  cpc  ,  one  obtains, 


V  •  Ax  =  ■ 


-q  v  (x-vt) 


4718  c 


(Hy/cY 


3/2 


(X_Vt)  ^  2^  2 
— - —  +  y  +z 

l-(v/c) 


.3/2  ’ 


(37) 


and 


d<Pc 

fit 


q  v  (x-vt) 


4718 


(i-(v/c)2 


3/2 


(X_Vt)2  2  2  ^ 
- — v+y  +z 

vHv/c)  j 


3/2  * 


(38) 


In  this  case,  JT  is  a  non-solenoidal  true  current  (because  it  is  due  to  a  single  moving 
charge)  and  fi(V  •  Ax)/fit=V*E[  ±  0  as  the  reader  can  verify  by  taking  the  time 
derivative  of  Eq.  (37). 

Comparing  Eqs.  (37)  and  (38)  shows 

V  •  Ax  +  (fi  cpc  /  fi  t)  /  c2  =  0  .  (39) 


Approved  for  public  release;  distribution  is  unlimited. 


8 


Eq.  (39)  closely  resembles  the  Lorenz  condition,  Eq.(12)  which  is  used  to  obtain  the 
primed  scalar  wave  equation  (Eq.  (18)).  There  are  important  differences,  however:  i) 
unlike  the  Lorenz  condition,  which  holds  for  the  net  primed  vector  potential,  A' ,  Eq.  (39) 
applies  to  only  one  component  of  the  net  vector  potential  A,  so  it  does  not  violate  Eq. 
(33);  ii)  Eq.  (39)  is  an  equation,  not  a  “condition”  to  be  met  by  adjusting  variables;  and, 
iii)  Eq.  (39)  is  in  covariant  fonn  in  the  unprimed  variables,  which  contrasts  with  the 
Lorenz  condition.  The  Lorenz  condition  does  not  apply  in  the  unprimed  variables,  a  fact 
that  is  generally  overlooked  in  comparisons  with  results  of  special  relativity. 

The  fact  that  an  equation  that  is  similar  to  the  Lorenz  condition  can  be  obtained  in  the 
unprimed  variables  is  well  known.  For  example,  Panofsky  and  Phillips  [1]  apply  the 
continuity  condition,  JT  +  dp  /  d  t  =  0  (which  refers  to  a  true  current  in  a  finite  length 

conductor,  with  a  time  varying  charge  density  at  the  end).  They  then  apply  Eq.  (36), 
giving  (without  subscripts), 


'[JT(r,  t')] 
r(t)-r  (t) 


dv 


(40) 


along  with  an  analogous  expression  for  the  scalar  potential, 


9c  (EO 


Af[pAlldv, 

4718  J  r(t)-r  (t) 


(41) 


to  show  that  Eq.  (39)  is  a  result  of  the  continuity  condition.  What  is  generally  overlooked 
it  that  the  vector  potential  in  this  case  is  the  component  related  to  the  true  current  rather 
than  the  net  vector  potential.  We  use  subscripts  to  preserve  the  distinctions,  which,  as  we 
have  stressed,  makes  all  the  difference. 


We  obtained  Eq.  (39)  without  invoking  the  continuity  condition;  however,  if  one  views 
current  in  a  wire  as  the  flow  of  a  line  of  charges  with  drift  velocity  v  over  a  stationary 
line  of  opposite  charges,  then  Eqs.  (34)  and  (35)  can  be  viewed  as  a  similar  effect  where 
the  dynamic  potentials  are  related  to  the  development  of  an  excess  charge  at  the  end  of 
the  conductor. 


A  final  point  is  that  some  of  the  more  subtle  features  of  relationships  among  component 
variables  tend  to  be  ignored  in  the  standard  approach.  Eq.  (39)  shows,  for  example,  that 
in  the  presence  of  time  varying  Coulomb  fields,  every  point  in  space  with  non-zero 
6  cpc  /  d  t  acts  as  a  point  source  contribution  to  AT  (see  Eq.  (29)).  While  this  component 

of  At  cannot  contribute  to  a  magnetic  field  B,  its  time  variation  does  contribute  an 
induced  field  E[  which,  in  principle,  can  be  detected  by  the  force  it  produces  on  charged 
particles. 
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The  Dynamic  Gauss’  Law 


In  the  standard  approach,  the  general  expression  for  the  divergence  of  the  combined 
electric  fields  in  the  presence  of  dynamic  Coulomb  fields  is  derived  as  shown  for  Eqs. 
(15)  and  (16).  All  approaches  simply  apply  the  static  form  of  Gauss’  law  (Maxwell’s 
Eq.(l)  )  to  the  sum  of  the  dynamic  variables,  Ec  and  E, ,  giving  the  expression, 

V«(Ec+Ej)  =  p/e.  (42) 

No  conceptual  basis  for  the  new  dynamic  variables  is  provided.  Despite  this,  Eq.  (40)  is 
generally  accepted  throughout  the  physics  literature,  presumably  because  it  provides  the 
standard  gauge  transformed  scalar  wave  equation  (Eq.  (18)).  Without  it,  one  cannot 
obtain  the  total  electric  field,  E.  We  now  show  that  Eq.  (42)  and  its  equivalent,  Eq.  (16) 
are  incorrect. 

First,  assume  for  the  moment  that  Eq.  (40)  is  valid.  The  original  V*E®  (Eq.(l))  and 
V  •  Ej  (Eq.  (22))  are  not  required  to  apply  in  Eq.  (42).  Instead,  the  divergence  of  the 

sum  of  the  two  fields  at  a  given  point  in  space  is  equal  to  the  charge  density  at  that  point. 
This  conflicts  with  the  Faraday  law:  the  complete  Faraday  law  must  always  apply,  so  Eq. 
(22),V*Ej=0,  always  applies.  Consequently,  Maxwell’s  Eq.  (1)  holds  for  both  the 

static  and  dynamic  Coulomb  fields  so  that  there  is  no  need  for  the  concept  of  dynamic 
fields  if  Eq.  (42)  is  correct.  The  source  of  the  difficulty  is  that  the  standard  assumption 
that  Maxwell’s  Eq.  (1)  applies  to  the  sum  of  Ec  and  Ej  is  simply  invalid. 

A  proper  derivation  of  the  dynamic  form  of  Gauss’  law  is  available  from  the  retarded 
field  expressions  for  cpc  and  AT  (Eqs.  (34)  and  (35)).  For  the  sake  of  brevity  we  do  not 
include  the  rather  lengthy  equations  that  result  from  the  straightforward  derivatives,  but 
give  only  a  summary  of  the  results,  which  the  reader  can  readily  verify.  If  one  obtains 
Ec  from  -  V(pc  and  E,  from  -SA,  /fit outside  the  singularity,  one  finds  that  V  •  Ec  ^  0 

andV  "E,1  ^  0  . 

On  the  other  hand,  if  one  now  uses  Eq.  (34)  and  (35)  to  compute  the  divergence  of  the 
sum  of  these  two  vectors  outside  the  singularity,  one  obtains, 

V  •  (Ec  +  E^)  =  0 .  (43) 

In  order  to  include  the  singularity  at  the  location  of  the  charge,  q,  at  (x-vt),  y,  z,  we  use 
the  fact  that  the  correction  term  yj  1  -  (v/c)2  is  only  valid  distances  larger  than  the  size  of 
the  charge  [3],  Retardation  effects  cannot  exist  in  the  immediate  vicinity  of  the  charge,  so 
the  fields  there  follow  the  static  inverse  square  law  for  E*  .  The  result  is  a  correct 
expression  for  the  dynamic  form  of  Gauss’  law, 
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V*(Ec  +  E*)  =  q/s. 


(44) 


Applying  the  principle  of  superposition  for  a  distribution  of  charges  as  in  the  static  case 
gives  the  more  general  expression, 

V«(Ec  +  Ej’)=p/e.  (45) 

This  result  reveals  a  previously  unrecognized  feature  of  the  dynamic  Gauss’  law,  namely, 
that  it  holds  for  the  dynamic  Coulomb  field  and  for  only  one  component  of  the  induced 
field,  E[ .  There  are  no  contributions  from  the  induced  fields  generated  by  the  two 
displacement  currents. 

Returning  to  the  topic  of  variable  labels,  note  that  Maxwell’s  Eq.  (2)  holds  for  general 
Coulomb  fields.  Thus,  it  applies  to  both  the  static  E*  and  to  the  dynamic  Ec  described 
here. 


The  Scalar  Wave  Equation 

The  correct  form  of  the  dynamic  Gauss’  law  reveals  a  flaw  in  the  gauge  approach.  The 
key  step  in  obtaining  the  scalar  wave  equation  relies  on  the  incorrect  form  of  Gauss’  law, 
(Eq.  (16)).  We  now  see  that  there  is  no  such  difficulty  with  the  unprimed  variables  and 
the  correct  version  of  the  dynamic  Gauss’  law. 


Rewriting  Eq.  (45), 


-V2(pc  -d(V •  At)/ <3t=p/s . 


(46) 


Applying  Eq.  (39)  to  Eq.  (46)  gives, 

V2cpc  - fi2cpc  / <3 12  / c2  =  -p/s  ,  (47) 


which  is  the  scalar  wave  equation  in  unprimed  variables.  Eq.  (47)  is  independent  of  the 
vector  wave  equation.  This  contrasts  with  the  results  of  the  gauge  approach,  where  the 
primed  scalar  and  vector  wave  equations  are  actually  meaningless  by  themselves.  (For 
example,  even  in  cases  where  actual  Coulomb  fields  are  absent,  the  gauge  approach 
requires  that  one  consider  an  imagined  Coulomb  field.)  Both  primed  equations  are 
coupled  because  only  the  sum  of  the  primed  electric  fields  is  accessible  to  laboratory 
testing.  A  further  advantage  of  Eq.  (47)  and  the  vector  wave  equation  (derived  in  the  next 
section)  over  that  obtained  from  the  gauge  approach  is  that  they  can  legitimately  be 
discussed  in  terms  of  the  unprimed  variables  that  can  be  measured  directly  or  compared 
with  results  from  special  relativity. 
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The  Vector  Wave  Equation 


We  begin  with  Eq.  (6),  which  is  Maxwell’s  Eq.  (5)  in  terms  of  the  vector  potential, 

Vx  Vx  A  =  V(V •  A)- V2A  =  p  Jt  +  ps<9Ec/ St  +  ps<3Ej/dt .  (6) 

Eq.  (6)  shows  there  are  three  source  tenns  for  the  net  vector  potential,  A.  To  underscore 
the  point  that  A  is  actually  the  sum  of  components,  we  consider  each  separately.  The 
component,  Ax  ,  due  to  the  true  current  has  already  been  discussed.  The  explicit  retarded 

field  expression  for  the  vector  potential  component  due  to  the  Coulomb  displacement 
current  can  be  written  as 


[sSEc(r',f)/ot] 

r(t)-r'(t) 


dv’. 


(48) 


and  the  third  component  of  A,  due  to  the  displacement  current  from  the  total  induced 
field,  Ej ,  can  be  written  as 


[8SEi(r',t')/St]dv, 

r(t)-r'(t) 


(49) 


The  net  vector  potential  is  the  sum  of  these  three  terms, 

A=AX  +Aj  +Ac  .  (50) 

Since  the  source  current  for  net  A  is  solenoidal,  the  general  criterion  for  S(V  •  A)/o  t  =  0 
is  satisfied.  Taking  the  time  derivative  of  both  sides  of  Eq.  (6),  and  employing 
d(V  •  A)/dt  =  0 ,  gives  the  wave  equation  for  the  induced  field,  E( , 

V2E,  -  psd2Et  /dt2=  p  <9Jt  /  5 1  +  ps  32Ec  /fit2.  (51) 

Thus,  with  the  present  approach,  one  obtains  exactly  the  same  wave  equation  for  E ,  as 

that  obtained  directly  (Eq.  (21))  from  the  basic  set  of  Maxwell’s  equations.  As  with  the 
scalar  equation,  the  vector  wave  equation  holds  for  the  unprimed  variables.  This 
consistency  check  is  not  possible  with  the  gauge  approach  since,  by  definition,  the 
meanings  of  the  primed  variables  differ  from  those  of  the  unprimed  variables. 

For  completeness,  we  note  that  the  corresponding  wave  equation  for  B  is  obtained  from 
the  curl  of  both  sides  of  Eq.  (6).  Using  Maxwell’s  Eq.  (2)  and  the  fact  that 
V  x  V(V  •  A)=0  gives, 
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v2B-nsd2B/at2=  nd(VxjT)/at, 


(52) 


which  is  identical  to  that  obtained  directly  from  Maxwell’s  Eqs.  (2),  (3),  and  (4). 

Another  consistency  check  is  that  the  solenoidal  source  terms  in  both  Eqs.  (51)  and  (52) 
generate  the  solenoidal  vector  fields,  B  and  Ej .  The  same  is  not  true  of  Eq.  (6)  since  both 
sides  of  Eq.  (6)  are  solenoidal,  regardless  of  V  •  A  .  On  that  topic,  we  have  shown  that  the 
basic  wave  equations  (i.e.,  Eqs.  (49),  (50),  and  the  scalar  wave  equation  (Eq.  (47)))  can 
be  derived  from  the  potentials  without  invoking  an  electromagnetic  gauge.  The  correct 
wave  equations  are  obtained  for  any  time  independent  function  for  V  •  A  .  The  reason  can 
be  seen  from  Eq.  (29).  A  constant  V  •  A  generates  neither  a  magnetic  field  nor  an  electric 
field,  so  it  is  physically  undetectable.  It  is  physically  meaningless.  Furthermore,  there  can 
be  no  physical  justification  for  the  existence  of  such  a  term  where  solenoidal  current 
sources  generate  continuous  flow  lines  for  A.  So,  the  concept  of  a  source  tenn  for  A  is 
physically  meaningless,  which  requires  V  •  A=  0  .  Equation  (6)  can  then  be  written  as, 

V2A-fi2A/5t2  =-pJT-psdEc/at.  (53) 

As  a  final  point,  the  present  results  reflect  a  general  underlying  physical  principle: 
solenoidal  current  sources  generate  solenoidal  fields  so  the  statement  that  there  are  no 
monopoles  for  B,  may  be  extended  to  include  E(  and  A. 


Reconciling  with  Special  Relativity 

The  direct  incorporation  of  gauged  (or  transformed)  variables  into  the  equations  for 
special  relativity  is  not  correct  because  one  cannot  equate  primed  and  unprimed  variables. 
We  already  noted  that  AT  =  v(pt ,  where  is  the  component  due  to  the  true  current,  and  that 
the  unprimed  analog  to  the  Lorenz  condition  is  also  expressed  in  terms  of  AT .  Bearing 
this  in  mind,  if  one  re-examines  the  primed  vector  potential  wave  equation  (Eq.  (14))  it  is 
clear  that  it  also  applies  only  to  AT .  The  displacement  contribution  from  the  Coulomb 

field  is  removed  in  the  gauge  transformation  process  leaving  a  non-solenoidal  source 
tenn  and  a  wave  equation  for  a  non-solenoidal  vector  field.  In  other  words,  the  primed 
net  vector  potential  is  the  same  as  the  component  due  to  the  true  current,  A'  =  AT  .  Thus, 
the  main  result  of  the  standard  gauge  in  electromagnetism  is  a  wave  equation  that  is 
applicable  to  only  AT  .  Since  the  standard  covariant  expressions  are  based  on  the  primed 

vector  fields,  it  is  clear  that  all  such  expressions  actually  refer  to  the  AT  component, 
rather  than  the  net  A.  Hence,  simply  changing  the  label  of  the  A  vectors  in  the  special 
relativity  equations  to  AT  corrects  the  misleading  equations  and  removes  the  conflict 

between  the  present  results  and  the  requirements  of  special  relativity.  Reiterating,  the 
present  results  show  that  Maxwell’s  equations  are  consistent  with  V  •  A=() ,  while  the 
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relativistic  equations  are  consistent  with  V  •  AT  +  (d  cpc  /  d  t)  /  c2  =  0  . 

A  similar  situation  exists  in  the  quantum  mechanical  treatment  of  electromagnetism. 
Feynman  and  Hibbs  [4],  for  example,  discuss  the  fact  that  a  quantum  mechanical 
formalism  that  incorporates  the  Lorenz  condition  is  considerably  more  cumbersome  than 
that  using  V*A=0.  So  V  •  A=0  is  generally  assumed  to  apply.  Again,  the  conflict 
between  this  requirement  and  special  relativity  is  resolved  by  recognizing  that  the 
relativistic  equations  refer  to  AT  rather  than  A. 


Summary 

We  have  shown  that  much  of  the  standard  formalism  of  electromagnetism  cannot  survive 
the  application  of  precise  definitions  to  its  variables.  The  key  to  addressing  this  problem 
is  the  inclusion  of  an  additional  equation  to  the  basic  set  of  Maxwell’s  equations.  This 
additional  equation  is  required  to  complete  the  statement  of  Faraday’s  law.  With  a 
complete  set  of  basic  equations,  all  the  essential  relationships  needed  for  valid  derivations 
of  the  wave  equations  can  be  obtained  without  an  electromagnetic  gauge. 
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Appendix 


In  the  present  paper  we  attribute  errors  in  the  standard  derivations  of  wave  equations  to 
the  common  practice  of  neglecting  distinctions  among  key  field  variables.  Reference  [5], 
describes  similar  errors  in  a  variety  of  textbook  analyses  related  to  applications  of  the 
Faraday  law.  This  Appendix  provides  a  simple  illustrative  example  of  such  errors  using  a 
circular  shell  where  no  Coulomb  fields  exist.  This  example  also  serves  to  support  several 
of  the  central  points  made  in  the  main  text. 

Consider  a  thin,  conducting  circular  cylindrical  shell  of  1  meter  diameter,  with  a  loop 
resistance  of  1  ohm,  enclosing  a  uniform  magnetic  field  which  is  increasing  linearly  in 
time  so  that  the  rate  of  change  of  flux  (emf)  is  1  volt  (Figure  1).  (The  shell  is  thin  enough 
that  skin  effect  can  be  neglected.)  What  is  the  potential  difference  between  points  A  and 
B? 

A  common  approach  is  to  apply  the  Faraday  law  in  the  following  manner:  Using  the  fact 
that  the  line  integral  around  any  closed  path  is  equal  to  the  negative  of  the  flux  enclosed, 
one  conventionally  assumes:  i)  that  no  fields  exist  inside  a  good  conductor  and  ii)  any 
convenient  path  can  be  selected  to  give  the  field  between  points  A  and  B.  The  simplest 
choice  is  a  closed  path  that  passes  along  a  zero  field  circumferential  segment  from  A  to  B 
inside  the  shell  and  exits  the  shell  at  B  and  follows  the  diameter  back  to  A  to  close  the 
loop.  This  loop  encompasses  half  the  area,  so  the  non-zero  field  between  A  and  B  is  0.5 
volts/meter  and  the  potential  between  A  and  B  is  0.5  volts.  If  one  considers  the  symmetry 
of  the  setup,  however,  this  must  be  incorrect. 

The  root  of  the  problem  is  ignoring  the  distinctions  between  Coulomb  fields  and  induced 
fields.  Potentials  are  obtained  from  line  integrals  of  Coulomb  fields  while  emf  s  are 
obtained  from  line  integrals  of  induced  fields.  Furthermore,  paths  cannot  be  chosen 
arbitrarily:  the  induced  fields,  E, ,  actually  form  closed  loops  that  are  concentric  with  the 
conducting  shell,  (and  inside  the  shell  itself).  Because  perfect  circular  symmetry  is 
assumed,  there  can  be  no  induced  charges,  so  there  are  no  Coulomb  fields;  hence,  there 
are  no  potential  differences.  So,  there  are  no  fields  directed  along  the  diameter  between 
points  A  and  B.  Furthermore,  the  fields  inside  the  shell  are  not  zero,  so  a  1  amp  current 
will  be  generated.  The  final  point  is  that  the  physics  of  the  problem  cannot  described  by 
the  sum  of  Et  and  Ec  as  required  in  the  gauge  approach. 

This  example  also  relates  to  requirement  for  a  full  characterization  of  Ej  and  the  possible 
presence  of  a  source  term  for  E:  in  the  absence  of  a  Coulomb  field.  A  source  tenn  implies 
a  potential,  cp,  (See  Eqs.  (24)  through  (28)).  Thus,  the  line  integral  of  E,  between  any  two 
points  on  the  circular  shell  would  include  a  contribution  from  the  gradient  of  (pt .  The 
result  would  be  an  accumulation  of  charges  at  different  points  along  the  shell,  which,  in 
turn,  would  produce  detectable  potential  differences.  (A  concrete  example,  consistent 
with  the  circular  symmetry  of  the  present  example,  is  a  source  tenn  for  E;  along  the 
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central  axis  of  the  shell.  This  would  produce  a  potential  difference  between  the  inner  and 
outer  surfaces  of  the  shell.)  Such  sources  produce  no  effect  on  the  total  current,  however, 
since  the  line  integral  of  a  gradient  o f  cp ,  around  a  closed  circuit  is  zero.  The  presence  of  a 

source  term  would  yield  deviations  between  predicted  and  measured  potentials  in  a  wide 
range  of  electromagnetic  devices.  In  this  case,  one  can  fairly  assume  that  the  absence  of 
evidence  of  a  source  tenn  is  evidence  of  its  absence,  so  Eq.  (22)  must  apply  for  closed 
current  loops. 


A 


Figure  1:  Conducting  shell  enclosing  a  uniform,  time-varying  magnetic  field. 
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